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The statistics of critical spin-spin correlation functions in Ising systems with non-frustrated dis- 
order are investigated on a strip geometry, via numerical transfer-matrix techniques. Conformal 
invariance concepts are used, in order to test for logarithmic corrections to pure power-law decay 
against distance. Fits of our data to conformal-invariance expressions, specific to logarithmic correc- 
tions to correlations on strips, give results with the correct sign, for the moments of order n = — 4 
of the correlation-function distribution. We find an interval of disorder strength along which correc- 
tions to pure-system behavior can be decomposed into the product of a known n-dependent factor 
and an approximately n-independent one, in accordance with predictions. A phenomenological fit- 
ting procedure is proposed, which takes partial account of subdominant terms of correlation-function 
decay on strips. In the low-disorder limit, it gives results in fairly good agreement with theoretical 
predictions, provided that an additional assumption is made. 
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I. INTRODUCTION 

Magnetic systems with quenched disorder exhibit 
properties which often differ from those of their homo- 
geneous counterparts. The Harris criterion jjj for the 
relevance of disorder, with respect to possible changes of 
universality class at the magnetic phase transition, boils 
down to the question of whether the pure-system specific 
heat exponent a is positive or negative. In the former 
case, the universality class changes, relative to that of 
the pure system, while in the latter it does not. For 
the two-dimensional Ising model, the Harris criterion is 
indecisive, as the pure-system specific heat diverges log- 
arithmically. Thus, alternative schemes must be formu- 
lated in order to investigate this borderline case. Nowa- 
days, mounting evidence 0, 0| suggests that the critical 
behavior of two-dimensional Ising ferromagnets is only 
slightly modified by the introduction of non-frustrated 
disorder. The corresponding changes are given by loga- 
rithmic corrections to pure-system power-law singulari- 
ties. Though the theoretical framework underlying such 
corrections is well-understood, their detection in numer- 
ical studies turns out to be a subtle matter 0.0.0.01 

Upon introduction of disorder, the scaling dimensions 
associated to the moments of correlation functions at crit- 
icality acquire multifractal behavior Q , characterized by 
a set of independent exponents, as opposed to the pure 
case where a single exponent is expected. 
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It was analytically predicted by Ludwig @ that the 
moments of the spin-spin correlation functions of the dis- 
ordered Ising model behave asymptotically as: 



i?-"/ 4 (lni?) 



n(n-l)/8 



(1) 



where n = 1,2,..., and the label Jy denotes average over 
quenched disorder (bond) configurations. The power-law 
term i?~™/ 4 corresponds to the uniform-system correla- 
tion decay, while the logarithmic corrections are given by 
(lni?)"^ n 1 '' s , Critical correlation functions are non-self 
averaging quantities, that is, the width of their distribu- 
tion stays essentially constant as the number of samples 
N increases; however, the sample-to-sample fluctuations 
of averaged values (e.g., the assorted moments just men- 
tioned) do go down approximately with the square root 
of sample size as the latter increases 

For the two-dimensional Ising model, the "typical", or 
most probable, correlation function, defined for a specific 
disorder configuration (i.e., the zeroth order moment of 
the distribution), is predicted 0] to behave as: 



G = cxp [InG (i?)] av ~ R- 1 ' 4 (lni?) 



-1/8 



(2) 



where, again, the uniform-system behavior (G(R) ~ 
R~ v , r\ = 1/4) is accompanied by a logarithmic correc- 
tion. Note that corrections are expected to be absent for 
n = 1, i.e., for the average correlation decay. For higher 
moments n > 2, the correction terms should be present, 
with corresponding exponents A n = n(n — l)/8, given in 
Eq. GJ. 

Here we apply numerical transfer-matrix methods, to- 
gether with finite-size scaling and conformal invariance 
concepts, to investigate the random-bond Ising model on 



2 



strips of a square lattice (L x TV, N — > oo). Our purpose 
is to confirm the analytical results derived by Ludwig, 
for the exponents r\ n — n/A and A n = n(n — l)/8, asso- 
ciated to the moments of critical correlation functions of 
the model. 



II. METHOD 

We consider a random-bond Ising system on a square 
lattice, whose Hamiltonian is given by 



7/ = ~X] J v u i a i 



(3) 



where <Ji = ±1 are the (site) spin variables, and Jy 
are ferromagnetic interactions between nearest-neighbor 
spins, extracted from the quenched probability distribu- 
tion: 

P(Jij) = \ Wa - J o) + KJij - rJo)] , < r < 1. (4) 

The exact critical temperature j3 c = l/fcsT c is known 
from duality [HG3 



sinh(2/3 c J ) sinh(2/3 c r J ) = 1. 



(5) 



In our calculations we initially consider a system with r = 
1/4, for which Eq. © gives T c (l/4)/J = 1.239- •• (in 
these units, the critical point of the uniform system is at 
T c (l)/Jo = 2.269 • • • ). We shall also investigate stronger 
disorder, by using r = 0.1, 0.05, and 0.01 (for the lat- 
ter, the critical temperature is T c (0.01)/Jo = 0.5089 • • • ). 
The calculation of spin-spin correlation functions follows 
the lines of Section 1.4 of Ref . with standard adapta- 
tions for an inhomogeneous system^4[- Taking two spins 
on, say, row 1, separated by a distance R, and for a given 
configuration C of bonds, one has: 



(6) 

where <jq = {<Tq . . . aft} and correspondingly for ctr; the 
bonds that make the transfer matrices % belong to C. 
For pure systems the 2 L — component vectors ip, ip are 
determined by the boundary conditions; for example, the 
choice of dominant left and right eigenvectors gives the 
correlation function in an infinite system^]. Here, one 
needs only be concerned with avoiding start-up effects, 
since there is no convergence of iterated vectors. This is 
done by discarding the first few hundred iterates of the 
initial vector vo. From then on, one can shift the dummy 
origin of Eq. 10 along the strip, taking ip (ip) to be the 
left- (right-) iterate of Vq (vo ) up to the shifted ori- 
gin. In order to avoid spurious correlations between the 
dynamical variables, the respective iterations of -0 and ip 
must use distinct realizations of the bond distribution. 



As regards conformal invariance, this has been a fun- 
damental tool to the understanding of the critical behav- 
ior of pure systems in two dimensions. Specifically, the 
spin-spin correlation function on a strip can be obtained 
from the conformal mapping w = x + iy = (L / 2tt) hi z , 
of the plane z = u + iv onto a strip of width L with peri- 
odic boundary conditions across, leading to the following 
asymptotic behavior at criticality |l,'il |: 



^»purc 



tt/L 



(sinh 2 (7rx/i) + sin 2 (7ry/L)) 



1/2 



(7) 



where 77 = 1/4 for pure Ising spins. 

For disordered spin systems at criticality, one expects 
conformal invariance to be preserved when averages over 
disorder are taken. Indeed, numerical evidence for this 
has been found in several instances, e.g., random-bond 
Isi ng m odel [4L HoL ll4T | . random-bond g-state Potts mod- 
els random transverse Ising chains at T = (equiva- 
lent to the two-dimensional McCoy- Wu model) [13], and 
Ising spin glasses [H El S E3 El ■ 

We first recall the simpler case where correlation de- 
cay in the bulk disordered system is purely polynomial. 
Considering correlation functions on a strip, preservation 
of conformal invariance in the sense just defined means 
that the decay-of-correlation exponents rji , associated to 
the assorted moments rm of the probability distribution 
(PDF) of correlation functions in the disordered magnet, 
can be determined via an extension of Eq. (0, i.e. by 
assuming: 



m, r-j z m , z — (sinh 2 (7ra;/L) + sin 2 (7ry/i)) 



1/2 



(8) 



The rji can be extracted by fitting numerical data to 
Eq. iJHJ. Thus far, it has been found that results derived 
in this way (e.g. for spin glasses) are consistent with 
those coming from different analytical methods which do 
not explicitly invoke conformal invariance [Isl. I2IL 1221 ] . 

Before going further, one must recall that an extensive 
analysis of the inaccuracies associated to the finite strip 
width L, carried out in Ref. l2lL indicates that finite- width 
effects are essentially subsumed in the explicit L depen- 
dence given in Eq. Q. Thus, higher-order finite-size cor- 
rections presumably do not play a significant role in the 
conformal-invariance properties of assorted moments of 
the correlation-function PDF (at least when critical cor- 
relation decay in the bulk is purely power law-like). 

Turning now to the present case, one must envisage 
how the multiplicative logarithmic corrections, predicted 
in Eqs. and (j2j for bulk correlationSjjtranslate onto a 
strip geometry. According to Refs. US, 13, for distances 
R ^> L along the strip the correlation decay should be 
exponential, with m n = G n (R) ~ exp(— R/£ n ), where 

-4- =Vn-K7rbg(lnL) + 0(g 2 ) . (9) 

In the above, & is a universal amplitude related to the 
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normalization of the three-point function 12.1 124 
g(lnL) = 



and 



.9o 



1 + 7r b go In L 



(10) 



where go is proportional to the intensity of disorder. In 
the limit L — ► oo, Eq. turns into a go-independent 
form, 



Art 
hl~L 



InL' 



(11) 



A semi-quantitative verification of Eq. 11 111 was given, for 
n = only, in Ref . H 

In Sec. IIIII we at first apply Eq. © to our data, i.e. 
logarithmic corrections are not explicitly considered. The 
idea is to keep as close as possible to the framework in 
which conformal invariance is known to apply to disor- 
dered systems, and check whether the resulting effective 
exponents, rff^ , deviate from their pure-system values 
in a way consistent with the corrections predicted in 
Eqs. (QJ and. {pj. Then we attempt to extract quantita- 
tive information, by comparing our numerical results to 
Eq. ©, as well as its asymptotic form, Eq. 11 111 . Finally, 
we investigate a phenomenological scheme to incorporate 
subdominant corrections to exponential correlation decay 
on a strip, which is subsequently applied to the analysis 
of low-disorder data. 



III. NUMERICAL RESULTS 

Estimation of uncertainties associated to the several 
moments of the correlation-function PDFs follows the 
procedure described in Ref. flfl where it was shown that 
for finite- width strips of length N, the width of the PDF 
remains constant as TV varies. Though the dispersions 
AG (or A(lnG)) do not shrink down to zero, it is possi- 
ble to extract reliable information about average values, 
as the dispersion of these among independent samples 
is proportional to 1/y/N. One can check, for instance, 
that fluctuations of said averages are < 1% for a strip of 
length N — 10 6 . 

In Fig. we present the moments m n of the 

correlation-function PDF against the variable z, where 
we used L = 16 and N = 10 5 for 1 < x < 5, < y < 8. 
From least-squares fits of the data, one can find estimates 
for the effective exponents r]f s in the form suggested by 
Eq. . Results are displayed in Table [I] Deviations 
from pure-system values are consistent with the signs pre- 
dicted for logarithmic corrections from Eqs. (QJ and J2J, 
namely an enhancement for n — 0, approximate neutral- 
ity for n = 1, and reduction for n > 2. 

In order to produce a more quantitative analysis, we 
now resort to Eqs. (©, CU- 

In all cases to be described below, we used L = 10, and 
the horizontal distance x was taken long enough in order 
to try and detect the effects which give rise to logarith- 
mic corrections (2^,01) see Eqs. 10, 111 111 . In practice, 
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Figure 1: (Color online) Double-logarithmic plot of mo- 
ments of the correlation-function distribution, against z = 
(sinh 2 (7rx/L) +sin 2 (7ry/L)) 1/2 , for r = 1/4, L = 16, N = 
10 . Lines are least-squares fits of single power-law form to 
data, from whose slopes one extracts estimates for the rj^ of 
Eq. ©. 



Table I: Estimates of effective exponents rji , from least- 
squares fits of moments of correlation-function distributions 
to the form rm ~ z Vi . Data for r — 1/4, L — 16 and 
N = 10 5 . 



0.271(2) 
0.252(2) 
0.473(4) 
0.671(8) 
0.85(1) 



we used 12 < x < 35 lattice spacings along the strip, 
as it was noticed that the use of larger x gave rise to 
large sample-to-sample fluctuations, thus compromising 
the accuracy of our data. Also, and for the sake of sim- 
plicity, we made y = 0. 

Here we evaluated 10 independent estimates for N — 
10 5 for the moments of the correlation-function PDF, 
and took error bars as three times the standard devia- 
tion among these. This was necessary, since it is these 
uncertainties which constitute the single major impedi- 
ment to the accurate estimation of the exponents rj n and 

An- 

In the semi- logarithmic plot of Fig. one can see that 



4 




i i i I i i i i I i i i i I i i i i I i 



10 15 20 25 30 
R 



Figure 2: (Color online) Semi-logarithmic plot of moments 
of the correlation-function distribution, against distance R 
along the strip (i.e. making y — in Eq. 0). Straight lines 
are least-squares fits of data, from whose slopes the assorted 
values of A Table HH have been estimated (see text). Here 
r = 1/4, L = 10, N = 10 5 . 



Table II: Estimates of exponents A„, from least-squares fits 
of moments of correlation-function distributions to a pure- 
exponential decay form, as given by Eq. Hill , keeping the ??„ 
fixed at their pure-system values: 770 = 1/4, rj„ — n/4 (n > 1). 
Data for r = 1/4, L — 10 and N = 10 5 . Uncertainties in last 
quoted digits are shown in parentheses. 



n 


A n (calc.) 


An (expected) 





-0.05(1) 


-1/8 


1 


-0.0012(7) 





2 


0.081(1) 


1/4 


3 


0.218(2) 


3/4 


4 


0.405(4) 


3/2 



an exponential decay with distance, which is a precondi- 
tion for Eqs. (JJ, 111 111 to apply, has set in to a good degree 
of accuracy for the range of R depicted. We first checked 
whether the values of strip width L = 10, and disorder 
r = 1/4 used so far, might correspond to the asymp- 
totic regime described by Eq. (11 II . The corresponding 
estimates for A„, shown in Table ITT1 have been extracted 
from least-squares fits of our data to the form given in 
Eq. Ill II . by keeping the rj n fixed at their pure-system 
values: r/o = 1/4, r) n = n/A (n > 1), 

As can be seen, though the estimates extracted from 
our data have their sign as predicted, their magnitude 



is at most only about one-third of the expected (except 
for n = 1, in which case the calculated result may be 
taken as broadly consistent with the predicted absence 
of corrections) . The most likely explanation for such dis- 
crepancy is, therefore, that the asymptotic regime has 
not been reached yet. Similar considerations were used 
in Ref. 23, where the 4-state Potts model was analysed 
with help of the generic (i. e., non-asymptotic) form, 
Eq. <PJ. From free-energy data, the numerical values of 
g{\nL) for L = 6 — 10 were extracted, and it was shown 
that they explained the respective deviations of gap scal- 
ing data from their known exact values. 

A major obstacle to following the same steps here is 
the presence of randomness, which severely limits the 
accuracy of both free-energy (see, e g., Ref. Q) and 
correlation-function estimates. Furthermore, the exact 
bulk free energy is not known, contrary to the case of the 
4-state Potts model 0. 

We then compared correlation-function data pertain- 
ing to different values of L, namely L = 10 and 16 (both 
for r = 1/4). We found that, e. g., fitting L = 16 results 
to the asymptotic form Eq. l(TTTl gives estimates for A n 
which fall essentially within the error bars of those given 
in Table [H] for L = 10. In addition to the randomness 
effects just mentioned, such poor resolution also reflects 
the logarithmic L-dependence of the corrections under 
study, and indicates that for the values of L < 25 within 
current reach of transfer-matrix calculations, there is lit- 
tle hope of extracting reliable information from varying 
L. 

Next, we investigated the behavior of our results 
against varying disorder. This means varying r in Eq. 
which corresponds to changing go in Eq. 11 Oil . The actual 
functional dependence go = go(r) is not known, apart 
from the basic fact that go must increase as r decreases 
towards zero. A similar degree of freedom is absent in 
the 4-state Potts model investigated in Ref. US From 
Eqs. ©, dl , and itTTll . by doing this for fixed L one 
expects to move along the crossover between small- and 
large-L regimes. 

We calculated moments of order 0—4 of the correlation- 
function PDFs, for r = 0.1, 0.05, and 0.01. Fits of our 
data to the asymptotic form, Eq. ltTT)l . are shown in 
Fig- El For each n^l, the sign of calculated exponents 
remains in agreement with predictions and their abso- 
lute value increases with increasing disorder (Ai remains 
close to zero, reaching 0.07(1) for r — 0.01). However, 
the trend against increasing disorder is clearly towards 
overshooting the theoretical predictions, with no sign of 
stabilization. 

It seems reasonable to exclude the possibility of non- 
monotonic behavior for even smaller r, which would then 
bring the estimated exponents back to their predicted 
values. The simplest explanation for the behavior de- 
picted in Fig. is that, at r < 0.05 or thereabouts (for 
fixed L = 10), the effects of the percolation fixed point at 
r = 0, T c (0) = 0, are beginning to manifest themselves. 
Indeed, the description synthesized in Eq. Q implicitly 
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Figure 3: (Color online) Logarithmic correction exponents 
A n , calculated via the asymptotic expression, Eq. Qllfl from 
moments of the correlation-function distribution, against dis- 
order intensity r (disorder increases with decreasing r, see 
Eq. (Jlji). Error bars are of the same order of, or smaller than, 
symbol sizes. Dashed horizontal lines starting on vertical axis 
mark expected values of X„ = —1/8, 0, 1/4, 3/4, and 3/2 
respectively for n = - 4. All for L = 10, TV" = 10 5 . 



envisages a perturbation scenario where pure-Ising crit- 
ical behavior is modified by a marginally-relevant oper- 
ator, which induces the logarithmic corrections investi- 
gated here. Thus it is plausible that such framework 
would fail to include the extreme disorder close to the 
zero-temperature fixed point. 

In order to check this hypothesis, we must go back to 
the non-asymptotic form, Eq. QJ, and ask whether (for 
fixed L) there is a range of disorder in which the correc- 
tions to the pure-system exponents can be decomposed 
into the product of an n-dependent and an n-independent 
factor (respectively, A„ and nbgQ.nL)) as in Eq. ©, and 
whether such range is limited by increasing disorder. As 
pointed out above, in contrast with Ref. [23 here we are 
allowed to vary go, which can then substitute for the L- 
variation investigated in that Reference. 

We assumed the r\ n and A n to have their exact pre- 
dicted values and, from the numerical estimates of mo- 
ments of correlation-function PDFs, calculated 7r6g(ln L) 
for n — 0, 2, 3, and 4. The results are exhibited in 
Fig. 31 We first note that the n = data show a much 
wider range of variation than the rest. This may be re- 
lated to specific properties of the zeroth moment of the 
correlation-function PDF, in whose calculation the oc- 
currence of small values plays a disproportionately im- 
portant role (compared to the moments of order n > 1). 



Figure 4: (Color online) Estimates of irbg(lnL) of Eq. IjllH 
(assuming rj n , A n to have their exact values) calculated from 
moments of the correlation-function distribution, against dis- 
order intensity r. Error bars are of the same order of, or 
smaller than, symbol sizes. All for L = 10, N = 10°. 



However, we cannot advance a specific explanation at 
this point. We decided to concentrate on the analysis of 
n = 2 — 4 data, whose average gives itbglhiL) — 0.30(3) 
for r — 1/4, and 0.56(8) for r = 0.1. For smaller r we 
see that the scatter becomes larger, in a way consistent 
with the idea that the physical picture of Eq. © becomes 
frayed as disorder increases. The agreement is qualitative 
and, to a fair extent, quantitative, inasmuch as both Fig- 
ures and concur in indicating r ss 0.1 as the threshold 
beyond which the perturbative view definitely fails (for 
L m 10, that is; for larger L such threshold would pre- 
sumably move towards higher disorder). 

We now return to the relatively low-disorder case of 
r = 1/4, and investigate whether a perturbative scheme 
exists, which may account for the features exhibited by 
correlation functions on a strip geometry. Specifically, 
Eq. © predicts that, in the disordered system, corre- 
lations on a strip still decay purely exponentially, as in 
pure-Ising strips, only with a small change in the expo- 
nent. In the crossover regime to which L = 10, r = 1/4 
certainly belong, such incipient modification may turn up 
as an effective subdominant (i.e. less than exponential) 
correction to our numerics. Thus, one could expect an 
effective behavior: 

G n {R) ~cxp(-^R?7„/L)[i?] A " , (12) 

where g n are the pure-system exponents, and subdomi- 
nant corrrections are assumed to take the simplest pos- 
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sible form of a power law, with the A„ to be deter- 
mined. Following this reasoning, and recalling the be- 
havior of bulk correlations invoked in Eqs. (QJ and (2J, 
one sees that a literal transposition of the -R-dependence 
of the latter onto the z-dependence of strip correlations 
(in analogy with Eq. JHJO, would mean: 

G n (z) ~ [lnz] A " = [z (In z) 9 "p" , (13) 

where rjo = 1/4, A = -1/8, r) n = ra/4, A„ = n(n - 
l)/8 (n > 1), and q n = — A n /?7 n . In the regime x = 
R L, y — considered in Eq. O, which corresponds 
to z ~ exp(7ri?/L) in Eq. (§J, this indeed amounts to 
writing 

G„(i?)~exp(-7ri?r/ n /L)[i?] A " , (14) 

i.e., pure-Ising exponential decay with power-law correc- 
tions. Since Eq. (T3j) does not, as far as we are aware, 
have a rigorous justification on the basis of conformal in- 
variance, the predicted identification A„ = A„ extracted 
from comparison of Eq. fl2} to Eq. (Hj) must be regarded 
warily. Nevertheless, from the preceding analysis summa- 
rized in Tables U and [HI for each n one expects A n and 
A„ to have the same sign and, possibly, to be of similar 
orders of magnitude. Therefore, in what follows we shall 
take considerations based on Eq. ((13)) as a starting point 
to reanalyze the set of data displayed in Fig. [21 

Using Eq. Ill -it , the moment of order n = must corre- 
spond to q = 1/2. In Fig. fitting data for the moment 
of order n — we at first fixed q = 1/2 (as suggested by 
Eq. (2J), and allowed 770 to vary in Eq. 11 -111 . The best 
fit was found for 770 = 0.2505(7), corresponding to a x 2 
per degree of freedom (Xd f ) equal to 0.012. We then 
fixed 770 at the central estimate just found, and allowed 
Ao of Eq. i[T2l to vary. The Xd o f f° r this latter fit are 
displayed in the insert of Fig. with a clear minimum 
at A = -0.125. 

According to Eq. the moment of order n = 1 is 
not expected to present logarithmic corrections. In order 
to check this, we fixed q — in Eq. ltl3ll and allowed 
771 to vary. The best fit, shown in Fig. OH was found 
for 771 = 0.2505(1), very close to 1/4 of a uniform sys- 
tem, corresponding to a \d o f = 0.02. In the insert of 
the Figure, one can see that, upon keeping 771 fixed at 
the central estimate quoted above, and allowing Ai of 
Eq. 1121 to vary, the closest fit indeed occurs for Ai = 0. 

When turning to higher-order moments n > 2, we no- 
ticed that it was not possible, in general, to produce good 
fits of the whole range of data displayed in Fig. [21 to the 
variables suggested by Eq. ill .311 . However, subsets of the 
data did fit well against such variables, as is illustrated 
in Fig. We shall return to this point at the end of the 
Section. 

For the moment of order n = 2, again we fixed 
q = —1/2 in Eq. (13)) and allowed r\i to vary. The 
best fit was found for r) 2 = 0.4992(8), corresponding to 
a Xd o f = 0-03. Data for this setup are displayed in the 
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Figure 5: (Color online) Main diagram: double-logarithmic 
plot of moment of order n = of correlation-function 
PDF, against [2 (biz) 9 ] for q = 1/2. Line corresponds to 
770 = 0.2505 (best estimate obtained for fixed q = 1/2, see 
text); L = 10, r = 1/4. Insert: semi-logarithmic plot of 
Xd.o.f. against — Ao, for fits of n — moment data against 
z-' va [in z] A °, with 770 = 0.2505 (fixed). 



Table III: Estimates of exponents r) n and A„, from least- 
squares fits of moments of correlation-function distributions 
obtained with G„(z) ~ [z (In z) q "]~ Vrl , as given by Eq. 11311 : 
see also Eq. (TJJ. Data for r = 1/4, L = 10 and N = 10 s . 
Uncertainties in last quoted digits are shown in parentheses. 



n 


7?n 


A n 





0.2505(7) 


-0.124(4) 


1 


0.2505(1) 


0.000(2) 


2 


0.4992(8) 


0.248(6) 


3 


0.751(2) 


0.75(1) 


4 


1.032(6) 


1.54(4) 



main diagram of Fig. 03 In the insert of the Figure, we 
kept 772 fixed at the central estimate just found, and al- 
lowed A2 of Eq. ll"12ll to vary. The Xd o f °f the respective 
fits are shown, with a clear minimum at A2 = 0.25. 

The previous procedure was applied for the moments 
of orders n — 3 and n = 4, The results for all investigated 
exponents, r\ n and A„ are displayed in Table ITTTI 

It can be seen from Table ITTTI that both accuracy and 
precision (the latter, upon comparison with theoretical 
predictions of A„) of results become lower for higher- 
order moments; in particular, our estimate of 774 al- 
ready differs somewhat from 774 = 1 predicted in Eq. (21 
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Figure 6: (Color online) Main diagram: double-logarithmic 
plot of moment of order n = 1 of correlation-function PDF, 
against [z (lnz) 9 ] for q = 0. Line corresponds to r/i = 0.2505 
(best estimate obtained for fixed q = 0, see text); L — 10, 
r = 1/4. Insert: semi-logarithmic plot of Xd.o.f. against Ai 
, for fits of n = 1 moment data against z~ vi [lnz] Al , with 
Vi = 0.2505 (fixed). 




z(ln z)" 1/2 



Figure 8: (Color online) Main diagram: double-logarithmic 
plot of moment of order n = 2 of correlation-function PDF, 
against [z (lnz) 9 ] for q = —1/2. Line corresponds to 772 = 
0.4992 (best estimate obtained for fixed q = —1/2, see text); 
L — 10, r — 1/4. Insert: semi-logarithmic plot of Xd.o.f. 
against A2 , for fits of n = 2 moment data against z~ V2 [In z\ A2 , 
with 772 = 0.4992 (fixed). 
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Figure 7: (Color online) Double-logarithmic plot of higher- 
order moments of correlation-function PDF, against [z (In z) q ] 
with q = —1/2, —1, and —3/2 respectively for n = 2, 3, and 
4. Thick continuous lines show range of data used in least- 
squares fits of the rjn for q fixed as above. L = 10, r = 1/4. 



(though that for A4 is consistent, at the margin, with 
A4 = 3/2 from the same Equation). For n>4we have 
found that the situation worsens considerably. This is 
because: (i) short-distance deviations from asymptotic 
behavior tend to increase with n, and (ii) at the rel- 
atively long distances where asymptotic behavior sets 
in, the correlation functions themselves, and a fortiori 
their higher moments, become smaller, and their aver- 
ages prone to higher relative fluctuations. We recall the 
trend illustrated in Figure against increasing n: for 
all three sets of data shown there, the thick continuous 
lines (corresponding respectively to 19 < x < 32 [n = 2), 
21 < x < 32 (n = 3), and 20 < x < 31 (n = 4)) de- 
note the range of data used to extract the results quoted 
in Table HHI One sees that, since \q\ increases with n, 
sets of data for the same range of x become more com- 
pressed when plotted against [z(lnz) 9 ], thus increasing 
uncertainties upon fitting. 

As one might expect, such a picture of good fits to 
Eq: (1121 . producing estimates of r\ n and A n quite close 
to the exact values, does not hold for higher disorder. 
For instance, already for r = 0.1, n = 2, our best fit to 
Eq. O gives 772 = 0.433(1), A 2 = 0.019(5) { X 2 dof = 
0.04). Fixing 77 2 = 1/2, A 2 = 1/4 gives xlo.f. « 24 - 
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IV. CONCLUSIONS 

We have tested the theoretical predictions, given in 
Eqs. (QJ and (3J, of logarithmic corrections to the 
asymptotic behavior of spin-spin correlation functions, in 
the two-dimensional (unfrustrated) random-bond Ising 
model. The use of a strip geometry has been ac- 
companied by the corresponding finite-size scaling and 
conformal-invariance considerations. So far the validity 
of the latter, when dealing with non-uniform systems, 
refers to disorder-averaged quantities, and has been es- 
tablished only in cases where critical correlation decay is 
purely polynomial fiH l2ll l22T| . 

From examination of short-distance data (see Fig. 
and TableQJ , we could only establish that the effective ex- 
ponents obtained were consistent with corrections of the 
same sign as predicted by theory. The situation is similar 
to that found in Ref. 0, where the conformal mapping of 
the square lattice on a finite square was used, and only 
n = was investigated . Those authors concluded that, 
for the size of lattices used, they would be able only to 
assert that the sign of corrections was as predicted. 

By considering suitably long distances along the strip, 
we attempted to reach the regime where the effect of log- 
arithmic corrections on a strip geometry is expected to be 
more easily isolated, as predicted in Refs .l2li.l24L Allow- 
ing the intensity of disorder to vary produced a picture 
in which one is able to identify the app roximate range of 
validity of the results given in Refs . 123.1241 even for strips 
of relatively small width. This relies essentially on iden- 
tifying the interval of disorder strength r, along which an 
approximately n-independent factor can be isolated from 



numerical data, in accordance with the non-asymptotic 
forms predicted in Eq. jsj. The corresponding data are 
shown in Fig. Q] 

We also produced a phenomenological fitting scheme, 
expected to work at low disorder, and synthesized in 
Eq. 11211 . from which estimates of the postulated expo- 
nents A„ were extracted. While the analogy to Eq. JHJ, 
used to establish Eqs. (|l3l) and ifplj) . seems short of a 
rigorous justification on the basis of conformal invari- 
ance, the ensuing prediction that A„ = A„ appears to be 
backed to a large extent by our numerical data for low 
disorder r = 1/4. As expected, the validity of such a 
scheme does not extend far towards increasing disorder. 

In summary, the evidence shown here provides a check 
of the validity of Eqs. (QJ and Q for disordered Ising sys- 
tems, in the appropriate limit of long distances. It would 
be interesting to ascertain whether any deeper connec- 
tion lies between the phenomenological form Eq. l(T2|l and 
conformal-invariance properties. 
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